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ABSTRACT 


The Fast Fourier Transform (FFT) has dramatically changed 
the Signal processing field by allowing the transformation 
of data from the time domain to the frequency domain at 
speeds rivaling and in many cases exceeding real time. The 
output spectrum obtained from the FFT depends on the ampli- 
tudes and the phase coherence of the sSinusoids that comprise 
the input waveform. It is the purpose of this investigation 
to modify the FFT so that it will rely more heavily on the 
phase coherence of these sinusoids rather than on their am- 
plitude. In this way a second indicator of signal presence 
or absence would be available to add more reliability to the 
estimation process. The study proposes a method for reduc- 
ing the dependence of the output coefficients on input ampli- 
tudes and investigates its merits with several types of data. 
A method of displaying the phase coherence indicators is 


presented. 
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I. INTRODUCTION 


Signal processing deals with the problem of detecting 
Signals in the presence of noise. If the noise is station- 
ary GausSian with a mean value of 0, and if the signal is 
absolutely stable, then it is theoretically possible to dif- 
ferentiate between the signal and the noise, no matter how 
small the signal to noise ratio. By integrating long enough 
(coherent integration) or averaging enough transforms (inco- 
herent integration) the noise will eventually average out to 
zero or achieve a state of zero variance. This is the abso- 
lute ideal and is never found in practice. The noise is 
rarely, if ever, purely Gaussian, nor is it absolutely sta- 
tionary. Signals are not absolutely stable, nor do they 
exist for infinite periods of time, thereby disallowing in- 
finite integration or averaging times. Because of these 
real-world constraints there are limitations to the effec- 
tiveness of all signal processing methods when applied to 
practical problems. 

The Fast Fourier Transform (FFT) is one of the most com- 
mon and effective methods by which data can be transformed 
from the time domain to the frequency domain. This trans- 
formation produces a frequency spectrum that gives Die rela- 
tive power in the various frequencies that comprise the time 
domain data. Due to the practical limitations of signal 
processing stated above this method of signal processing is 


limited in effectiveness. It is the purpose of this study 





to try to develop another transformation, similar to that of 
the FFT, that will produce a frequency spectrum that is more 
dependent on the phase coherence of the different frequency 
components that constitute the time domain data and less de- 
pendent on their power content. The spectrum analyst will 
then have a second indicator to aid him in making his deter- 
mination of signal presence or absence. Hopefully this will 
result in more reliable estimates and allow signal detection 
in some situations that have previously been beyond the phys- 
ical capabilities of the FFT. 

Part 1 of this investigation developes a procedure for 
obtaining a numerical indicator of phase coherence from the 
FFT. The procedure involves weighting the amplitudes of the 
components that add to produce each of the FFT output coeffi- 
cients equally so the resulting spectrum will rely more 
heavily on the phase relationships between additive components 
than on their amplitudes. Part 2 compares the merits of the 
procedure's signal enhancement capabilities with that of the 
FFT. Part 3 summarizes these comparisons and Part 4 draws 
some conclusions about the differences between the Modified 


FFT's spectrum and that of the FFT. 
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II. PART 1 DEVELOPMENT 


A. GENERAL 

The Fast Fourier Transform (FFT) is an algorithm for com- 
puting the Fourier Coefficients of a discretely sampled wave- 
form X(t ) where the subscript T implies X(t) to be finite in 
length (of length T) and also periodic in t with period T. 
f,=1/T is the fundamental frequence of the FFT. 

Appendix A gives a brief description of the Discrete and 
Fast Fourier Transforms. The flow graph of a 16 point FFT 
along with some of the notation to be used throughout this 
paper is presented in Figure 1-1. The term "level" will be 
used to denote the resultant vector formed after each pass 
(matrix multiplication) of the algorithm. Level (0) refer- 
ences the sampled version of the input waveform in its natu- 
ral order. level (1) references the input vector after bit 
reversal has taken place. Level (2) is the intermediate re- 
sult (or vector) resulting from the first pass (matrix multi- 
plication) through the FFT algorithm. The number 2 is used 
to imply that two elements of the input vector are involved 
in the computation of each element of the intermediate vector 
at level (2). The notation is similar for the remainder of 
Pnewoalcoritam. A Circle with a + or —- in it implies addition 
SieaWeuraction 1S to take place. A circle with a 1,2,3,... 
in it implies multiplication by W', W*, W°,... where 


W = exp(-j27™/N). 
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The magnitudes of the output coefficients from the FFT 
depend on the amplitude of the input signal and its phase co- 
herence throughout the transformation. When the FFT algo- 
rithm is viewed as a Sequence of vector (or phasor) rotations 
and additions, it is the phase coherence of a signal that 
determines how the intermediate vectors will add in produc- 
ing each output coefficient of the FFT. 

Each pass through the algorithm consists of a sequence of 
vector rotations, followed by a vector addition. As previous- 
ly stated, the length or magnitude of each output coefficient 
is a function of the magnitude of the intermediate vectors 
and of the manner in which they add to make up each output 
coefficient. The term phase coherence is used to describe 
the net effect of the phase relationships throughout the en- 
tire transform. A pure Sinusoid would have phase coherence 
throughout the entire transform, whereas the phase relation- 
Ship for noise would be completely random and incoherent. 

It is phase coherence that allows the intermediate vectors 
of a pure sinusoid to add end-to-end in a straight line, 
whereas the intermediate vectors for noise add at random 
phase angles, thereby producing output coefficients with an 
amplitude less than the maximum value they could achieve if 
they were to add ina straight line. 

To give a better understanding of the concept of phase 
coherence consider the vector additions in Figure 1-2. The 
resultant vector may be viewed as the output vector of the 


FFT and its component vectors as the intermediate vectors of 
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Vector addition analogy of FFT 


(a) Incoherent vector additions 

(b) Coherent vector additions 

(c) Incoherent vector additions, same as (a) 
but with intermediate components normalized 


to same length. 


Figure 1-2 
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the FFT that have been rotated and added to form the result- 
ant or output vector. Figure 1-2a could be viewed as inco- 
herent noise with no particular order in the manner in which 
the component vectors add to produce the resultant vector. 
Conversely, Figure 1-2b shows a resultant vector of what may 
be interpreted as the result of the transformation of a pure 
Sinusoid with zero .phase or frequency shift throughout the 
time record being transformed. Each component adds to the 
others in a straight line and at the phase angle of the input 
Sinusoid. The vector in 2a achieves its magnitude not 
through phase coherence but simply because one of its com- 
ponent vectors dominates the others in magnitude. If each 

of its component vectors were restricted to having unit 
length while still maintaining its same phase angle, then 

the magnitude of the resultant vector would be forced to de- 
pend on phase coherence (the ability of the component vectors 
to add in a straight line) throughout the transform and not 
on the magnitude of any of its component vectors. This sit- 
uation is depicted in Figure 1-2c. 

When the output of the discrete correlation performed by 
the FFT is squared, the resulting spectrum is proportional 
to the, power contained in the original waveform at each ot 
phew correlating frequencies n/T (n=0,1,2,-,-2,N/2-1) for real 
valued input waveforms and is referred to as the power spec- 
tral density of the input waveform. The association of the 
output coefficients of the FFT is strongly related to, and 


dependent on, the power of the input waveform at its various 
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frequencies, or more explicitly, at the correlation frequen- 
cies of the FFT algorithm. 

It would seem plausible that if the magnitude of the in- 
termediate vectors could be removed then the output spectrum 
would rely more on the sinusoidal nature of the input wave- 
form than on the power of the input sinusoid. Or conversely, 
the ability of large vectors, due to incoherent noise, to 
corrupt or mask a coherent sinusoid would be reduced. The 
objective is to develop an indicator of the presence of a 
Signal through phase coherence with as little emphasis as 
possible on Signal or noise power, an indicator which could 
be used in conjunction with the FFT to rerun estimating sig- 


nal presence or absence. 


B. PREMISE 

The basic premise of this study is to modify the FFT 
algorithm in such a manner that the output coefficients give 
‘an indication of the relative phase coherence of the sinu- 
soids that make up the input waveform. To test the modified 
transform to see if it meets the premise, sinusoids that 
complete an integral number of cycles across the time record 
will be considered, that is, sinusoids that are integral 
multiples of the fundamental frequency of the transform fo. 
These are of course the correlating sinusoids of the FFT algo- 
rithm. It will be required to show that the phase angle of 
the output coefficient corresponding to a pure sinusoidal 


morte at any frequency nf,(n=1,2,...,N/2-1) reflects the 
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phase angle of the input sinusoid. This must be true if the 
intermediate vectors are to add in a straight line for a 
phase coherent sinusoid. 

If the resulting output coefficient representing an in- 
put sinusoid of frequency nf, reflects the phase angle of an 
input sinusoid at these frequencies, then the premise has 
been satisfied. Sinusoids at frequencies other than integral 
multiples of f» are seen by the correlation process as a sum- 
mation of sinusoids at the frequencies nf). That is, the 
FFT output does not reflect the phase angle of the input 
Sinusoids at frequencies other than integral multiples of fy 
“but instead approximates these sinusoids by sinusoids that 
are at integral multiples of f,). Therefore the phase angles 
obtained do not reflect the phase angle of the input sinu- 
soid, but instead reflect the phase angles of the correlating 
Sinusoids that, when summed, will approximate the input wave- 
form. Therefore it is meaningful to test the premiSe on 


input sinusoids that are integral multiples of fo. 


C. IMPLEMENTATION OF THE MODIFICATION 

The output coefficients of the FFT are in general com- 
plex numbers of the form c=at+jb with magnitude |c| = (a2+b2)? 
and phase angle 0=Tan ‘(b/a). When an output coefficient of 
the FFT is viewed with the objective of determining the phase 
Poele of its corresponding input sinusoid, the phase angle 


obtained from 0=Tan ‘(b/a) is the average difference in phase 


angle between the input sinusoid and a reference cosine. If 


BET: 





a Sine wave is used for the reference sinusoid then the phase 
angle 6=Tan~ (a/b) must be used. In addition, one must know 
whether the correlating sinusoids were of the form cos(2mnk/N) 
~jsin(2mnk/N) or cos(2mnk/N)+jsin(2mnk/N). Either form re- 
sults in the same magnitude but with a different phase angle. 
When 6 is computed as Tan” !(b/a), changing the sign of the 
Sine wave correlation results in a reflection of the phase 
angle across the real axis. When 6 fe computed as Tan ?/(a/b), 
the sign change results ina reflection across the imaginary 
axis. The phase angle has physical meaning if one recognizes 
what the reference Sinusoid was and whether exp(-j2mnk/N) or 
exp(j2mnk/N) was used for the correlation process. 

After each pass of the FFT algorithm (at each level), the 
intermediate coefficients can likewise be represented in the 
same manner as the output coefficients. These intermediate 
coefficients are analogous to the intermediate vectors re- 
ferred to in Figure 1-2. The object will be to make these 
intermediate coefficients unit vectors, vectors of magnitude 
1. In this way the phase relationships of these vectors will 
dominate the length of any given output coefficient with as 
little emphasis on the magnitude of these vectors as possible. 
The method utilized will be simply to divide each coeffi- 
cient by its magnitude. This process will be referred to as 
normalization Since it normalizes the magnitude of the coef- 
frements to a value of 1 and in no way alters the ratio of 
faeereal to imaginary part of each coefficient, that is, the 


Pages angle ofeach coefficient is retained. 
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The problem then becomes one of determining where the 
normalization should take place in the FFT algorithm. The 
following properties of the FFT butterfly will be used to 
make this determination for real valued input waveforms: 

(1) <All of the coefficients at the 2 level are real 
valued since each coefficient at this level is either the 
sum or difference of 2 of the original real valued input co- 
efficients. 

(2) Half of the coefficients at the 4 level are complex 
numbers, namely the odd numbered coefficients x (keise5e 0 
N-1) and these can. have any phase angle, the rest are real 
valued. 

(3) The odd numbered coefficients at the 4 level affect 
only the odd numbered output coefficients (those associated 
with ff) and its odd harmonics), and similarly, the even num- 
bered coefficients affect only the even numbered output co- 
efficients (those associated with the even harmonics of fo). 

(4) Above the 4 level an increasing number of coeffi- 
cients become complex, for example at the 8 level the complex 
coefficients are: X (k=1,2,3,5,6,7,9,10,11,...,N-1). At the 
16 level the complex coefficients are: x Gea Ao Oe . 
eee 1 2543,14,15,.+.,N-1L), etc. 

Mathie input coefficients are normalized, the only pos- 
eiole phase angles that they can have are 0 or 180 degrees 
Mere with a mapnitude of 1) since they ane all real numbers. 
With N input coefficients then there will only be pe possible 


Normalized input vectors (if the value zero is allowed then 
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there will be 3N possible input vectors). This results ina 
finite number of output phase angles and thus by contradic- 
tion prohibits the vectors from adding in such a manner as 
to reflect all possible input phase angles. 

The same argument is true for the 2 level where the co- 
efficients again are real numbers. In fact the number of 
possible output phase angles is less than when the coeffi- 
cients were normalized at the input. The number of differ- 
ent vectors is still oN at the 2 level but there is one less 
pass between the 2 level and the output than between the in- 
put and the output, thereby reducing the number of passes in 
which coefficients are combined. 

The problem of determining where to implement this nor- 
malization procedure is greatly simplified if one remembers 
that the purpose of this modification is to achieve vector 
additions of the intermediate coefficients so that these vec- 
tors will add in such a way that the resultant vector is 
oriented in the direction of the input waveform's phase. In 
general that phase will be measured with respect to a cosine 
waveform throughout this study. 

What the above suggests is that a coefficient should be 
complex before it is normalized in order not to restrict the 
total number of phase angles achievable by the output coeffi- 
cients. 

The 4 level is the first level where complex numbers 
arise. As previously stated half of the coefficients affect 


oaiyeche output coefficients related to -even barmonrves=oe!l i942 
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At this level the transform can be split up into two distinct 
and independent parts, one for the odd harmonics and one for 
the even harmonics. The real coefficients are the building 
blocks of the even harmonics, by the previous argument of a 
finite number of output phase angles it will not be possible 
to normalize these coefficients and obtain all possible 
phase angles at the output. 

The complex coefficients at this level are the interme- 
diate coefficients that make up the output coefficients re- 
lating to the odd harmonics of fy. If these coefficients 
are normalized the output coefficients corresponding to the 
odd harmonics of fo» will reflect the phase angle of any in- 
put sinusoid at the frequencies f),3f),5f,),...,(N/2-1)f 9. 
Demonstration: 

Consider the bit reversed form of a 16 point transform 
of some function x(t) and the FFT butterfly for these coef- 
ficients up to the 4 level in Figure 1-3. As the notation 
indicates each coefficient at the 4 level is made up of four 
Sliements o& the input vectem esEach gmeuping of four c@efifi-— 


cients at the 4 level can be rewritten as: 


X"(4j) = x(O+j) + x(N/2+5) + x(N/44+5) + x(3N/443) 
K*(4j+1) = x(O+j) - x(N/2+5) + JIx(N/4+5) - x(3N/44+5)] 
K*(4j+2) = x(O+j) + x(N/2+5) - [x(N/44+5) + x(3N/44j)] 


K°(4j+3) = x(O+j) - x(N/2+5j) - JIx(N/4+j) - x(3SN/44))) 


merc j=0,1,2,...,N/4-1. This is a valid generalization for 
any N. The coefficients in question are the complex coeffi- 
cients and only these coefficients will be considered. Since 
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(0) (1) (2) (4) 

x(0) x(0) x(0)+x(8) x(0)+x(8)+x(4)+x(12) 

py) ~—C (8) x(0)-x(8) x(0)-x(8)+j(x(h)-x(12) ) 

x(2) = x() 5) a ware x(0)+x(8)-(x(h)+x(12)) 

x(3) x(12) x(4)-x(12) x(0)-x(8)- j(x(h)-x(12) ) 

x(h) sue ) ecg x(2)+x(10)+x(6)+x(1h) 

x(5) -x(10 IS (2) x(10)+5(x(6)-x(1h)) 
2 oink tacarnte (x06) toc ( 1s) 


x(6) x(6) 
x(6)-x(1h) x(2)-x(10)-j(x(6)-x(1)) ) 


x(1)+x(9) x(1)+x(9)+x(5) +x(13) 
x(1)-x(9) x(1)-x(9)+j(x«(5)-x(13)) 
mS) eds x(1)+x(9)-(x(5)+x(13)) 


x(5)-x(13) x(1)-x(9)-j(x(5)-x(13)) 


56'S) ciel ull) Kt eee 1) Exe oe 5) 
x(3)-x(11) x3 )—-x (11) 4 j (x07) -x(15)) 
x(7)+x(15) pall 3) awa ipa (a eee e105) ))) 


x7) —-x (ip x(3)-x(11)-j(x(7)-x(15)) 


ay th 


x(7) x 
me) (1 
x(9)  =-x(9) 
x(10) x(5) 
meri) 6x13 
rae) =X (3) 
x(13) x(11 
x(1h) x(7) 
cay x(15) 


K | 


XN 


The FFT algoritm for a 16 point transform up to the 4 level 


Figure 1-3 
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normalization is simply the division of any coefficient by 
its magnitude, then the normalization factor (the magnitude 
of the complex coefficients) is: 
{Ix(j) - x(N/2+j)]? + [x(N/44j) - x(3N/44+5)1?)}? 
The input sinusoid can be generalized as: 
x(t) = Asin(2tnfot + 9) where n=1,3,5,...,N/2-1, f 9=1/T 
and 9 is arbitrary ~1<0<T. 
The sampled version of x(t), namely TS) (k=0 45,2752 o7Ne)) 
is: x(k)=Asin(21mnfokAT + 6) where AT is the time increment 
between samples. 
Now: T=NAT and fo9=1/T=1/(NAT), then: 
x(k)=Asin(21nk/N + @). 
Therefore the divisor of the coefficients can be written, 
{ [Asin(2mnj /N+0)-Asin(27n(N/2+j)/N+6)]? 
+[Asin(2mn(N/4+3)/N+0)—-Asin(2mn(3N/4+5)/N+0)]2}2 
letting C=2t1nj/N and B=2t, then the divisor can be rewritten: 
A{[sin(C+0)-sin(C+B/2+6)]” “- [sin(C+B/4+6)-sin(C+3B/4+0)]2}2 
Now, Sin(xty )=sin(x)cos(y)+sin(y )cos(x) 
Expanding the above, 
A{{sin(C)cos(6)+sin(8)cos(C)-sin(C)cos(B/2+9)-sin(B/2+8 )cos(C)]? 
+[sin(C)cos(B/4+6 )+sin(B/4+6)cos(C) 
-sin(C)cos(3B/4+0)-sin(3B/4+8)cos(C)17}* 
Now, cos(xty )=cos(x)cos(y)-sin(x)sin(y) 


Expanding once more, 
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A{{sin(C)cos(6)+sin(6)cos(C)-sin(C)cos(B/2)cos(6) 
+sin(C)sin(B/2)sin(6)-cos(C)sin(B/2)cos(6) 
-cos(C)sin(6)cos(B/2)]? + [sin(C)cos(B/4)cos(@) 
-sin(C)sin(B/4)sin(6)+cos(C)sin(B/4)cos(6) 
+cos(C)sin(6)cos(B/4)-sin(C)cos(3B/4)cos(6) 
+sin(C)sin(3B/4)sin(6)-cos(C)sin(3B/4)cos(6) 
-~cos(C)sin(6)cos(3B/4) 12} 
Now B=2mn (n=1,3,5,...,N/2-1) 

=27,67,107,147,... 
then B/2=1,31,51,77,... 
and sin(B/2)=0 

cos(B/2)=-1 

Any emsan/ 2) 5m 2/2, oe 
Pigeebeeias3p 8 t1)7/21tt 
cos(B/4)=0 
3B/4=3n/2,90/2,157/2,210/2,... 
lyn 25) 24 sees Gel wets ae 

cos(3B/4)=0 
Making the above substitutions: 
A{[sin(C)cos(6)+sin(6)cos(C)+sin(C)cos(6)+cos(C)sin(6)]? 
+ [#sin(C)sin(@)+cos(C)eos(6)#sin(C)sin(8)#cos(C)cos(6)]2}” 
=A{[2sin(C)cos(6)+2sin(6)cos(C)]*+[{+2sin(C)sin(6) 
+2cos(C)cos(6)]2}* = A{4sin?(C)cos?(¢)+8sin(C)cos(6)sin(6)cos(C) 
+4sin*(8)cos?(C)+4sin’(C)sin’(6)~8sin(C)sin(®)cos(C)cos(8) 
+4cos?(C)cos?(0)}? = 2A{sin2(C)[cos?(6)+sin2(6)] 


eee (aun2(0)+cos2(0))}* = 2A{sin2(C)+cos2(C)}2 = 2A 
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Or, the divisor at the 4 level for the complex coefficients 
is dependent only on the amplitude of the input sinusoid; at 
is independent of N,n,j, or 8. This is equivalent to divid- 
ing all of the coefficients by a constant for any given in- 
put sinusoid of the form x(t)=Asin(2mnf»)t+8) and in no way 
alters the phase angle of the output coefficients. Its net 
effect is that of scaling the input sinusoid such that it is 
replaced by: 

x(t )=0.5sin(2tnfot+98) for any amplitude A. 

It can likewise be shown that a similar normalization at 
the 8 level results in division of the coefficients by 4A, 
at the 16 level the divisor is 8A and in general the divisor 
will be (L/2)A where L is the level at which the normaliza- 
tion was performed. 

The second half of the transformation can not be normal- 
ized at the 4 level for reasons already stated. At the 8 
level the coefficients once more divide, that is, the trans- 
formation can be broken down into 4 independent transforma- 
tions, each part being related to one-fourth of the output 
COemil@icntsme The intermediate coefficwents 1, 590R03 we... 
make up their corresponding output coefficients, likewise 
3,7,11,15,... make up their corresponding output coefficients. 
These coefficients offer no problem since they were already 
eomolex at the 4 level. Coefficients 0,4,3,72,... are still 
real valued at the 8 level and do not lend themselves to 
Hemmealization at this point. Coefficients 2,6,10,14,... 


were not complex at the 4 level but are complex at the 8 
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level. <A similar proof to the one used for normalization of 
the intermediate coefficients of the odd harmonics of fp at 
the 4 level can be used to Sironmtnest normalizing the coeffi- 
cients 2,6,10,14,... at the 8 level does not alter the phase 
angle of the output coefficients corresponding to these fre- 
quencies for sinusoidal inputs that are integral multiples 
of fo. Similarly at the 16 level the intermediate coeffi- 
cients that make up the frequencies nfy) (n=4,12,20,28,...) 
become complex and can therefore be normalized without alter- 
ing the phase angle of the output coefficients corresponding 
to these frequencies. 

This process continues until one reaches the output of 
the FFT at which time all of the coefficients (except F(0), 
the D.C. term) have become complex and therefore normaliza- 
tion could have been performed. The idea of normalizing dif- 
ferent coefficients at different levels produces one serious 
drawback. Normalization fixes the maximum magnitude an out- 
put coefficient can achieve. This is due to the fact that 
every intermediate coefficient somewhere in the FFT algorithm 
has been replaced with a magnitude of 1. Knowing where the 
normalization occurred one then knows how many more passes 
will be performed before the algorithm is complete. Since 
each pass consists of rotating and then adding coefficients, 
the number of vector additions remaining per output coeffi- 
cient is then equal to the number of remaining passes. Since 
the coefficients have magnitude 1, they can achieve a magni- 


tude at the output no larger than QP where P is the number of 
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remaining passes to be performed. For instance, in a 32 
point FFT, if the intermediate coefficients that make up the 
odd harmonics of fo» are CoH ReA MAR the 4 level there will 
be three remaining passes yet to be performed. At the 8 
level their maximum magnitude is 2, at the 16 level it is 4 
and at the 32 level (the output in this case) their magnitude 
can be 8. The only way the value 8 can be achieved is if the 
input was in fact a pure sinusoid at a frequency nf) (n=1,3, 
5,..-.).- It is important to note that if one knows how many 
points (N) were ina transformation, and what frequency a 
particular output coefficient is related to, that is, if one 
knows N, then the maximum value that coefficient can achieve 
is also known. From this information one may be able to 

make a prediction as to the probability of the existence of 

a signal at the frequency. If one also considers the magni- 
tude of this coefficient from more than one transform, then 
his prediction becomes more reliable. This in itself could 
be a powerful tool in signal processing work since the out- 
put coefficients are independent of the amplitude of the in- 
put signal. 

The serious drawback occurs when one analyzes a spectrum 
where the intermediate coefficients have been normalized at 
different levels. It is a simple procedure to multiply each 
output coefficient by the number of the level at which it 
was normalized. In the previous example a 32 point transform 
was considered and it was shown that the maximum magnitude 


that an output coefficient, dependent on normalization of 





its intermediate coefficients at the 4 level, could achieve 
was 8. Then 8 times 4 = 32. Likewise an output coefficient 
dependent on normalization of its intermediate coefficients 
at the 8 level can achieve a magnitude of 4. Then 4 times 

8 = 32. In general, the magnitude of the output coefficients 
can be made a constant (namely N) for an input sinusoid of 
the form Asin(2mnf,t+6) (n=1,2,3,...,N/2-1) independent of 

A, n, and 8. And the phase angle 6 oan be obtained from the 
output coefficients. 

When the input waveform consists of signal plus noise, 
the multiplication process leads to erroneous results. To 
best explain this problem, consider the output coefficient 
which is dependent on normalization at the yth level. This 
will always be the frequency (N/4)f>» Since this frequence co- 
efficient does not become complex until after the last pass 
of the FFT algorithm has been performed. This coefficient 
will always have a magnitude of 1 after normalization and it 
will have a magnitude N after it has been multiplied by its 
normalization level, which is always N in this case. That 
is, the frequency (N/4)fo will always have the magnitude N 
regardless of whether there was signal or noise at that fre- 
quency. To a lesser degree this same problem extends back 
to the frequencies that were normalized at the N/2 level, 
Poeiieaiay (0/8) £0 and (3N/8)f,. In general the frequen- 
cies that are normalized late in the transform tend to domi- 
Hate the spectrum. This can be explained if one considers 


that the multiplication factor was simply an artificial means 
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of adding att deere the output coefficients so that all 

of the frequencies could achieve the same maximum value. 
Frequencies that have their components normalized early in 
the transform rely more heavily on phase coherence to achieve 
large amplitudes. The frequencies that have their components 
normalized late in the transform have less time to achieve 
large amplitude through phase coherence. 

One solution to the above problem is to compare coeffi- 
cients that are dependent upon normalization at a given lev- 
el independent of the remaining coefficients. This has the 
drawback that the reliability of the estimate of signal pres- 
ence is reduced since the comparison is made against fewer 
coefficients. In addition it is inconvenient and lacking in 
continuity since at best every other coefficient is compared 
leaving gaps between frequencies. 

A second solution was devised in which all of the coef- 
ficients that depend on normalization at a given level are 
summed. Their mean value is determined and then each coef- 
ficient is divided by this mean value. This has the effect 
of giving the coefficients dependent on each level of normal- 
ization a mean value of 1. The spectrum can then be recom- 
bined with the effect of giving the entire spectrum a mean 
value of 1. A quick example will illustrate this procedure. 

Consider a 16 point FFT with an input signal consisting 
of equal amplitude i hertz and 2 hertz sinusoids plus noise. 
Figure 1-4 shows the input sinusoids and the sinusoids plus 


femse. Fipfure 1-5 shows the FFT algorithm and Figure 1-6 
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FFT 


FFT Normalized to mean of 1.0 
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Comparison of the output coefficients of the FFT 
and Modified FFT before and after they have been 


given a mean value of 1 


Figure 1-7 
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Shows the modified FFT algorithm starting from the 4 level 
(the FFT and modified FFT algorithms are the same up through 
the 4 level). Figure 1-7 compares the output coefficients 
of the modified FFT and the FFT before and after the coeffi- 
cients have been given a mean value of 1. The FFT coeffi- 
cients have also been given a mean value of 1 for comparison 
purposes ae Shows the value of this procedure for com- 
paring sequences of numbers when there is no known relation- 
Ship between the two sequences. In the case of the FFT the 
mean value of the coefficients was computed using all of the 
coefficients as opposed to the manner in which it was per- 
formed in the case of the modified FFT. There is no loss in 
generality in doing this because dividing a Sequence by its 
mean value in no way alters the relative magnitudes of the 
coeiticients. ) 

It is interesting to note in the case of the modified 
FFT that the signal frequence fo (the 1 hertz sinusoid) ap- 
pears much larger than 2f, before the "mean 1" process was 
applied. After applying this process the magnitude of fo and 
2f>5 are very close, as expected. Grouping the output coeffi- 
cients of the modified FFT that depend on normalization at 
each level and then comparing them with the FFT coefficients 
het. seen thateseae signals of the modified FFI have greater 
magnitude than the signals of the FFT and the noise coeffi- 


cients are smaller for the modified FFT than for the FIT in 


every case except 4fo. 
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LOT Mod FFT FFT Mod FFT 


f 9 L437 <1. UST 2f 9 ; 2000 < Lye 
3f 9 3405 2 OZ 6f y 7789 > .4756 
Of 9 (ico . 7184 

Tf o 1.3455 > 1.1494 4fo 4897 < 1.0 


The FFT spectrum was given a mean of 1 by dividing each 
coefficient by the mean of all of the coefficients. In the 
case of the modified FFT each coefficient was divided by the 
mean of the coefficients that were normalized at its normal-. 
ization level. Each procedure resulted in a spectrum with a 
mean of 1 but by different methods. To alleviate any doubt 
that the different methods used to give each process a mean 
value of 1 were responsible for the apparent improved proccss 
gain of the modified FFT over the FFT in this example the 
following compares the processes where each output spectrum 
has been given a mean value of 1 by the method used for the 


modified FFT. 


FFT Mod FFT FFT Mod FFT 
fo mezoor <= lLso0sr 2 9 ioe << eirozes 
oig fool” > ~o20u 6f o 7867 > .4756 
Of 9 .6804 < . 7184 


1.1494 Af o Js. 8) aoe. O 


V 


7TL£o 1.1860 


The method of normalizing coefficients to a mean of 1 
seems to have its merits but still does not solve the prob- 
lem of giving each frequency the same dependence on phase co- 


herence as opposed to amplitude plus phase coherence in the 
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case of the standard FFT. Consider the frequency (N/4)fo. 

It is in fact the FFT coefficient an a magnitude of 1. 

The frequencies (N/8)f,) and (3N/8)f») do not have their ampli- 
tude influences removed until just prior to the last pass 
through the algorithm. On the other hand the frequencies 

nfo (n=1,3,5,...,N/2-1) have their amplitude influences re- 
moved from their intermediate coefficients early in the 


mreansLorm. 


D. METHOD OF IMPLEMENTATION 

The method that was finally chosen was that of making all 
frequencies look like odd harmonics of 1/T.’ Keeping T con- 
stant and using frequency translation was tried but did not 
solve the problem as the even harmonics of f9, wnen trans- 
lated by 1/T, still behaved like even harmonics of fo when 
normalization was performed at the 4 level. That is, they 
did not reflect all possible phase angles at the output of 
the transformation. 

The method that did work was to vary T, the length of the 
transform, so that all frequencies were made to look like odd 
harmonics of 1/T. This is best explained by example. Con- 
sider a time record T seconds long and a sinusoid of fre- 
quency 2/T that has been sampled N times in T seconds. 

The correlating sinusoids are: 

cos(2mnk/N)-jsin(2mnk/N). 
The only non-zero correlation occurs when n=2. Now, if the 


time record were divided in half, then there would be N/2 
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samples in time T/2. The correlating sinusoids would then 
be: cos[2tnk/(N/2))-jsin[21mmk/(N/2)] and the only non-zero 
correlation would occur when n=1. The net result is that 

the correlating sinusoids are identical but the transform 
length is cut in half and the 2f) hertz sinusoid looks like 

a sinusoid of frequency fo to the time record of length T/2. 
If two transformations were performed and the output coeffi- 
cients are summed in the complex domain the net resuit is 

the same as if one transformation was performed of length T. 
This iS simply because the correlating sinusoids were identi- 
cal in both cases, that is, every multiplication was identi- 
cal in each case as was the summation of the products of 
these multiplications. If a 4f»5 Sinusoid were made to look 
like an odd harmonic, then the input record would have to be 
divided into four equal parts and 4-N/4 point transformations 
would then have to be summed. 

When the normalization modification is then added to the 
algorithm essentially only one-half of the algorithm is per- 
formed, the half that deals with the odd harmonics of fo. 
But the algorithm must be performed considerably more times. 
fancieald N poants are transformed to,obtain the frequency 
coefficients relating to nfy (n=1,3,5,...,N-1). Then 2-N/2 
point transforms are performed to obtain the frequency coef- 
ficients relating to nfo (n=2,6,10,...,N-2). Then 4-N/4 
point transforms are performed to obtain the frequency coef- 
miements relating to nfo (n=4,12,20,...,N-4), etce., until 
N/4-4 point transforms are performed to obtain the frequency 
Poetficient relating to (N/4)fo. 
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For a 1024 point transform the time required to obtain 
all of the frequency coefficients is approximately three 
times that of a normal FFT if the transformations are per- 
formed one after the other. If parallel processing were em- 
ployed the processing time could be reduced to half of that 
required for the FFT. If the criterion of NLog.,N (for the 
number of complex operations required in a transformation) is 
used for comparing the relative time nenuieeed for the two 
transformations then, for 32 points it requires 32 log,(32) = 
32x5 = 160 for the FFT. For the modified FFT it takes 


16x(5+4+3+2) = 16x14 = 224. For a 1024 point transform the 


FFT requires 1024x110 10,240 where the modified FFT re- 
quires 512x(10+9+...+2) = 512x54 = 27,648. For 8192 points 
the FFT requires 106,496 and the modified FFT requires 
368,640. The numbers do not account for the extra bit re- 
versals, additional ordering of output coefficients, or the 
normalization procedure required in the modified FFT. 

To summarize the results of this transformation: The 
proposed modification to the FFT algorithm will be more de- 
pendent on the phase coherence of a waveform to be trans- 
formed and less dependent on the power in the waveform than 
the FFT. When a sinusoid of the form Asin(2mnf)tt+6) (n=1,2, 
3,...,N/2-1) is transformed the resulting output coefficients 
will be a constant independent of A, n, or 8. And 0 Can be 
recovered from the output coefficient corresponding to n. 


The performance of this algorithm has been thoroughly in- 


Mestigated for signals of the form in the preceding paragraph. 
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Due to the nonlinearity of the process it is difficult to 
predict how it will perform in general. The next section in- 
vestigates the performance of the algorithm against that of 

! 


the FFT under identical conditions for different types of 


data. 
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III. PART 2 EXPERIMENTAL RESULTS 


A. MEASUREMENTS 

In determining the value of the modification to the FFT 
various types of data were transformed and compared with the 
FFT transformations. Several measurements were employed to 
evaluate the relative merits of the FFT modification when 
compared to the FFT. The output spectrum consisted of the 
magnitudes of the output coefficients referred to as (F, }, 
m-O)b.c,..,N/2-1). <A coefficient that is related to the 
frequency of an input signal is referred to as F.- In gen- 
eral, the objective of signal processing is the identifica- 
tion of a coefficient that is signal related from one that 
is not signal related. Normally this is accomplished by 
visual observation of the magnitude of the output spectrums 
of the transformations. Large coefficients tend to indicate 
the presence of a Signal and just as important, the consist- 
ency of the amplitude of an output coefficient from trans- 
form to transform indicates the presence of a Signal. It is 
the combination of these two effects that gives reliability 
to the estimate of signal presence or absence. | 

What may be one of the best means of determining the 
presence or absence of a Signal is the analysis by an exper- 
ienced spectrum analyst. Amplitude and consistency of an 
output coefficient do not always portray all of the informa- 
micen available in a coefficvent. There seem to be other 


characteristics in the behavior of a Spectrum from transform 
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to transform that gives extra support to an analyst's abil- 
ity to determine if a signal iS present or absent. In gen- 
eral it is difficult to explain what these additional 
characteristics are and it seems that only after a great 
deal of experience do they figure in to an analyst's logic 
when he decides if a Signal is present or not. The extra 
Signal characteristics generally pertain to a specific type 
of signal in a particular environment with which the analyst 
is very familiar; that is, he cannot generalize his experi- 
ence with a specific class of signals to all types of sig- 
nals and environments. 

The primary methods utilized in determining the merits 
of the modified FFT were: 

1. Estimation of signal enhancement through visual ob- 
servation. 

2. Determining the difference in amplitude between the 
magnitude of the signal and the magnitude of the mean value 
of the noise: 

n 
TE mat eae 


where Wn is the mean value of the noise coefficients. 


gL 


pemerl 


3. Determining the variance Of the noise: 
n 


ae. : oe - 2 
= pee! ts Wid Fos 


o2 
N 


normally associating the smaller variance with the better 


processing gain (referred to as smoothing the noise). 
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4. Utilization of a threshold detector to determine 
with what probability Fo would exceed all of the other EF or 
all of the FP. within a neighborhood of Poo such as Foose 
(i1=0,1,2,...,50), where the neighborhood in this case is 25 
coefficients on either side of Fo: 

DS. Comparing: m 7 

S n 

207 
for the different spectrums. The numerator relates how much 
ereater Fo is than the average of the comparison noise coef- 
ficients and the denominator relates to the smoothness of 
ioe noise coefficients. The form of this equation is identi- 
cal to that normally used in determining input signal to 
noise ratios. 

Method 5 was used in almost every experiment (it was not 
used in some of the early experiments) along with one or 
more of the other methods. It was found that methods 4 and 
Do agreed quite closely in almost every case (there were a 
negligible number of exceptions and then the difference it- 
Piowas teezligible). “Agreed” in the context of the above 
Sentence means that if a transformation was better than 
another in one test it was also better in the other test. 
Method 5 is the easiest to compute and relates directly to 
the signal to noise ratio of the output spectrum. AS a re- 


Sult it 1s used most often in the following pages to make 


comparisons of the FFT and the modified FFT. 
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The problem of determining input signal to noise ratios 
across a particular noise bandwidth was alleviated by com- 
paring the output spectrum of the modified FFT to that of 
the FFT. In this way the performance of the modification 
was judged relative to the FFT; i.e., the FFT was used as 
the standard of comparison. 

It was originally hoped that the resulting spectrum from 
the modified FFT would differ significantly from that of the 
standard FFT. If that was true, and if the modified FFT had 
sufficient process gain then the combination of the two 
spectra could result in a spectrum with a lower om and then 
hopefully a greater process gain. 

Since all of the theoretical work presented so far deals 
with simple sinusoids at multiples of the fundamentai fre- 
quency of the transformation a large void is left between 
what happens in that case and in the case of an input wave- 
form consisting of multiple sinusoids at arbitrary frequen- 
cies combined with random noise, Gaussian or not, especially 
Since the process is non-linear. 

Mathematical or statistical analysis of the modified FFT 
algorithm is complicated at best, and since the process is 
non-linear all conclusions as to its value were developed 


through experimentation with different types of input data. 
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B. SINGLE SINUSOID PLUS GAUSSIAN NOISE 

The first data chosen for analysis consisted of a simple 
Sinusoid combined with additive Gaussian noise. Most noise 
generators are designed to generate Gaussian noise and its 
statistical properties are well Known making it the most con- 
venient and generally accepted type of data for comparison 
purposes. Due to the simplicity of a single sinusoidal sig- 
nal and the availability of Gree ian medise generators the 
first experiments were performed on this type of data. 

The signal chosen was a 25 hz Sinusoid with additive 
Gaussian noise with a 100 hz bandwidth. The signal plus 
noise was generated with a "Wavetek" model 132 VCG/Noise Gen- 
erator. The input waveform was lowpass filtered at 100 hz 
and sampled at 256 samples per second. The transformations 
were 2048 points in length (8 seconds in time) with a re- 
sulting resolution of 1/8 hz. The noise to signal ratios 
ranged from O to 24 db in 3 db increments, for 3.33 minutes 
or 25, 8-second transforms at each level. 

Figures 2-1 through 2-9 show the results of these trans- 
formations. The horizontal axis represents frequency, in- 
creasing from left to right. 150 points (output coefficients) 
are displayed. The center of the highlighted area represents 
the 25 hz frequency. The left edge is 15.625 hz and the 
right edge is 34.375 hz. The vertical axis represents the 
relative magnitude of the coefficients. The fore and aft 
Sats represents timeswith the most current information being 


displayed at the front (or bottom) and the oldest information 
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being displayed at the back (or top). Each trace in each 
figure represents the average of ten transforms or 80 seconds 
of data. Each trace contains 80% of the same information 
contained in the traces adjacent to it. A boxcar type aver- 
ager is used with five compartments. In this case two trans- 
forms were averaged in each compartment, then the five 
compartments were averaged. Each succeeding trace is formed 
by dumping che compartment with the oldest information (20% 
of the information in the boxcar) and refilling it with new 
information, averaging and then displaying. With this in 
mind then only traces that are separated by five other traces 
are independent in terms of time from one another. 

Table I gives the signal to noise output ratios computed 
for the average of 25 transforms at each db level. The noise 
Statistics were computed from the 25 coefficients on each 
side of the signal coefficient. The coefficients immediately 
to the left and right of the signal were not included in 
these computations. As stated earlier all S/N output ratios 
were computed by subtracting the mean of the noise coeffi- 
cients from the signal coefficient and dividing the square 
of that result by twice the variance of the noise coeffi- 
cients. 

The most obvious result of comparing the respective out- 
put spectra of the FFT and modified FFT is their strong simi- 
larities. Close observation of the figures reveals the 
almost identical behavior of the respective spectra as the 


S/N ratio is reduced. Also the relative magnitude of the 
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Signal coefficient with respect to the magnitude of the co- 
efficients immediately adjacent to the signal coefficient 
seems to be the same for both spectra. This effect adds 
another element of similarity, the effect of having the same 
"shape". The observed variance of the noise appears to be 
somewhat greater in the modified FFT spectrum. 

Table I readily confirms the observed Similarities in 
the spectra. The S/N outputs for the two spectra are very 
close. The FFT generally has better S/N ratios and this can 
be attributed to the smaller noise variance observed in the 
photographs. The FFT becomes Significantly better than the 
modified FFT for the -18, -21, and ~-24 db S/N ratios, this 
is again confirmed in the photographs. Comparison of the 
FFT spectrum squared with that of the FFT spectrum multipliea 
by the modified FFT spectrum shows no signal enhancement 


over what can be gained by simply squaring the FFT spectrum. 
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Modified FFT 
25 hz signal, S/N = 0 db 
10 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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Modified FFT 

25 hz signal, S/N = -3 db 

10 transformed averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-2 








Modified FFT 


25 hz signal, S/N = -6 db 
10 transforms averaged per trace, 


transform length 2048 
Displayed bandwidth 18.75 Nay 0.12 


5 hz resolution 


Figure 2-3 
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Modified FFT 
25 hz signal, S/N = -9 db 
10 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-4 
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Modified FFT 
25 hz signal, S/N = -12 db 


10 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-5 
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Modified FFT 
25 hz signal, S/N = -15 db 
10 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-6 


O24 





- " 
AVEDA Wi bade 
wo ley ay 


“Whig We ian A : 


MON “rite 





Frit 


Ze 


pecans 


m 


re: AOI Se 





Modified FFT 
25 hz signal, S/N = -18 db 
10 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-7 
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Modified FFT 
25 hz sagnal, S/N = -21 db 
10 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-8 
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Modified FFT 


25 hz signal, S/N = -24 db 
10 transforms averaged per trace, transform length 2048 


Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-9 


ste) 


S/N out (db) 


ps FFT Mod-FFT 9 FFTXFFT nn ep 
O 26.295 27.260 39.790 40.902 
=€ 29.734 30.106 42.409 42.414 
os 35.007 32.860 48.469 47.258 
-9 S4.s20 sala7s 472235 47.578 
-12 25.480 24.053 32.632 32.208 
als 21.383 20.342 26.925 26.656 
-18 21.311 18.646 26.124 25.006 
-21 Lil, ale 7.664 13.726 12.490 
24 9.877 2.422 11.457 8.506 

Table I 


Average S/N out (db) vs. S/N in (db) for the Average of 25 
Transforms in the case of the 25 hz signal plus Gaussian 
noise. S/N out computed as: 

ae 2 

(Fi-u,) 
Zon 
n 

where F. is the magnitude of the signal coefficient and uF 


and oO, are computed from 25 noise coefficients on either side 


of F.. 
Ss 
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C. MULTIPLE SINUSOIDS PLUS GAUSSIAN NOISE 

To further investigate the relative merits of the modi- 
fied FFT Gaussian noise was again used. In this case multi- 
ple sinusoidal signals were present allowing investigation 
into the possibility of destructive interference between the 
signals due to the normalization procedure. 

The data came from a test tape made by the Naval Air 
Development Center. The resulting spectra are presented in 
Figures 2-10 through 2-14. The data consisted of four sig- 
nals in two groups of two. Onlyone group of two was consid- 
ered. The tape had ten minutes of data for each signal 
starting with three minutes of white noise followed by ten 
minutes of signals at -20 db and -18 db presented in close 
proximity to each other. These were followed by two more: 
Signals for ten minutes at -16 db and -14 db etc., until two 
Signals were presented at 0O db and +2 db (the last two were 
not transformed). The input signal to noise ratios were 
measured with respect to a 1 hz noise bandwidth. The fig- 
ures indicate where the signals are or where they should be. 

These figures are interpreted in exactly the same manner 
as Figures 2-1 through 2-9. In this case forty transforms 
of 2048 points each were averaged in each trace. The boxcar 
average was again used. This data was also filtered at 100 
hz and sampled at 256 samples per second. The frequency 
resolution is 1/8 hz. 

It should be noted that there was some question as to 


the reliability of the -18 db signal and it was not considcred 
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to be experimentally acceptable in the data analysis. The 
indication of where it should be in Figures 2-14 is only a 
best guess. , 

Table II gives the average computed output S/N ratios | 
for the summation of 40 transforms. These were computed 
using the same method employed throughout this study, In 
this case, when a Signal was on the right only 25 coeffi- 
cients to the right of that Signal were considered in com- 
puting the noise statistics. Likewise, if the signal was on 
the left only 25 coefficients to the left of that signal 
were considered in computing the noise statistics. In all 
cases the nearest coefficient to the signal coefficient was 
not considered in the computations. 

Table III is similar to Table II.except that 80 trans- 
forms were averaged instead of 40. 

The most important result of this experiment was to show 
that multiple signals did not interfere with each other and 
that the normalization process did not appear to destroy or 
reduce signal resolutions. Comparison of the figures again 
reveals similar spectra for the respective transformations. 
Comments identical to those made in Section III-B with re- 
spect to the single sinusoid plus Gaussian noise case could 
be made here with one exception. The spectra for the -16 
and -20 db signal to noise ratios tend to favor the modified 
FFT. In Figures 2-13 the top trace of the modified FFT 
shows the -16 db signal to be nearly as large as the -14 db 


Signal and in the second trace from the top it is actually 








larger. Also, the variance of the adjacent noise coeffi- 
cients appears to be smaller in the modified FFT spectrum 
than in the FFT spectrum. The individual distinctions are 
not as pronounced in the case of the -20 db Signal in Fig- 
ures 2-14 but the general appearance tends to favor the mod- 
ified FFT spectrum. 

Table II indicates that the FFT generally outperformed 
the modified FFT when the Signal to Rode ratios were rela- 
tively small. The reverse was true for the -16 and -20 db 
S/N ratios (as indicated above the -18 db signal is not being 
considered in the analysis). As in the case of the single 
Sinusoid presented in Section III-B no obvious trends were 
present in this data. Again the product of the two spectra 
when compared against the FFT spectrum squared did not prove 
to have any partilenier advantage. 

Table III compares the measured results when 80 trans- 
-forms were averaged instead of 40 as in Table II. The im- 
provement in S/N output was relatively consistent for all S/N 
ratios and for both spectra with neither spectrum gaining any 
particular advantage over the other as a result of the deeper 


average. 
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hO transforms averaged per trace, transform length 208 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
-2 db signal on right, -4 db signal on left 


peottwe = 282-93 UNIT oz INCH. 
fel @e = 3980-93 UNGTS7JNCH. 


MSOC Ti St TAPE UF TD 


94N= 72, 7-4 DB Figure 2-10 (FFT) 
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~l) db signal on left 
61 
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Figure 2-10 (Mod-FFT) 
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hO transforms averaged per trace, transform length 208 


Displayed bandwidth 18.75 hz, 0.125 hz resolution 


-2 db signal on right 
eee = JGR -4i UNIT 97 INCH. 


meets 


BAC TEST TAPE (MOD FFT) 


ei--2,-4 Db 
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4O transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
-6 db signal on left, -8 db signal on right 


Bee sLe = 39% -34 UNIT S7INCH. 
Meeohace = Je -23 LUNs TO INGH. 


NADC TEST TAPE (FFT) 


e7N=-5,-8 DB Figure 2-11 (FFT) 
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4O transforms averaged per trace, transform length 20,8 
Pasplayed bandwidth 18.75 hz, 0.125 hz resolution 
-~6 db signal on left, ~&8 db signal on right 


eee = Jee a UNITS INCH. 
Beeste — J.998-45 UNITS INCH. 


BoC IES! TAPE (MOD FFT) 


5 N= -§, -2 [BR Figure 2-11 (Hod-FFT) 
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4O transforms averaged per trace, transform length 20)8 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
mre ae Signal on left, -l2 db signal on right 


Beoodec = J.996-9) UNIT SZINGN, 
eorson = J. - 91 UNIT OMen. 


BNC iS] TAPE CFF?) 
o“N=~-10,-12 DB Figure 2-12 (FFT) 
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40 transforms averaged per trace, transform length 2018 
Displayed bandwidth 18.75 hz, 0. 125 hz resolution 
-10 db signal on left, <-1e ab Signal on right 


Beoealc = J04Gb-40 UNIT S7INCH, 
meocele = J2Q6C-48 UNGT o-JNCH. 


NADC TEST TAPE (MOD FFT). 
9/N= ~ 18, ae | Figure 2-12 (Mod-FFT) 
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hO transforms averaged per trace, transform length 208 
Displayed bandwidth 18.75 hz, 0. 125 hz resolution 
-il db signal on left, -16 ab signal on right 


Beet se = 3466-45 UNITS7 INCH. 
i= SCE 2 .SNGE -41 UNI TS7INCH. 


BC TEST TAPE (Fr!) 
N= - ie > UBD Figure 2-13 (FFT) 
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hO transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
-1h db signal on left, -16 db signal on right 


Peocalc = Je -93 UNITS7INCH. 
age = J895Ee oe) WNITOZINCH, 


BOC TEST TAPE (MOD FFT) 
S“N=-J4,-16 DB Figure 2-13 (Mod-FFT) 
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hO transforms averaged per Lae » transform length 2048 


Displayed bandwidth 18.75 hz, 0.125 hz resolution 
P10 db si@wal won right, ~20 ab signal on left 


Peockec ™ J35r-41 UNIT SZINCH, 
Beta ~ JM. -25 UNIT SINC. 


MC IEST TAPE (FFT) 
BN=- 16, =) NB Figure 2-ih (FFT) 
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hO transforms averaged per ate, transform ce 20,8 
Displayed bandwidth 18.75 hz, 125 hz resoluti 
18° db signal on right, -20 2 signal on left 


Meese = J.4992-93 UNITS/INCH. 
meets ct = 3.942 -4)] UNITS7INCH. 


BOC TEST TAPE (MOD FFT) 
BN-- 18. -OW DS Pigure 2-14 (Mod-FFT) 
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S/N out (db) 


a FFT Mod-FFT oo FFTxFFT PE 7% 
-2 19.655 18.059 23.985 23.165 
-4 17.086 16.388 20.218 20.092 
-6 2.064 5.086 0.358 2.521 
-8 5.901 8.397 4.622 7.018 

-10 6.882 5.693 4.769 5.373 

-12 8.215 7.121 6.761 6.641 

-14 8.847 7.258 7.499 7.160 

iG 5.489 6.386 4.008 ' 5.015 

-18 5.399 8.119 3.086 5.394 

-20 4.370 6.872 3.546 5.865 
Table II 


Average S/N out (db) vs. S/N in (db) for the Average of 40 
Transforms from the NADC Test Tape. S/N computed as: 
a 2 
(Fo-uL) 
Bo 
n 
where Ee. 1s the magnitude of the signal coefficient and Hy 


and oe are computed from 25 noise coefficients on one side 


ot F , 
Ss 
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S/N out (db) 


—— FFT Mod-FFT oo FFTxFFT etx 
-2 23.669 22.083 28.316 27.606 
' -4 . 19.943 19.738 23.094 23.300 
-6 3.334 ee ane 27 3.991 
-8 10.492 10.054 8.871 9.427 
= 8.575 8.366 6.420 7.140 
-12 10.395 8.889 8.827 8.437 
sil 10.895 10.048 9.414 9.021 
256 7.355 9.487 5.707 oe pl 
-18 7.504 11.319 4.676 7.294 
-20 7.032 10.372 6.065 8.581 

Table III 


Average S/N out (db) vs. S/N in (db) for the Average of 80 
Transforms from the NADC Test Tape. S/N computed as: 
7 2 
Ca ee, 
Boe 
n 
where ee is the magnitude of the signal coefficient and Hy 


and oO, are computed from 25 noise coefficients on one side 


oa F , 
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D. SINGLE SINUSOID PLUS GAUSSIAN NOISE PLUS BURST NOISE 

The data to be investigated in this section is similar 
to that discussed in Section III-B, except in this case addi- 
tional noise was added in the form of random bursts spaced 
at random intervals with random duration. The data was gen- 
erated in this manner to accentuate the effect of spectrum 
dominance by vectors with large amplitudes. This type of 
data, where large bursts of noise could dominate or corrupt 
the resulting spectrum, provides a good test of the FFT mod- 
ification. If the modification is performing as anticipated 
it will reduce the spectrum distortions caused by the large 
amplitudes of the burst noise by weighting all of the inter- 
mediate coefficients equally. Then large output coefficients 
will only be obtained through phase coherent vector additions 
and not through vectors with large amplitudes. 

In this case the signal was set at 10 hz vice 25 hz and 
additional noise was added in the form of bursts with random 
length and spaced at random intervals. Figures 2-15 through 
2-17 give samples of the input waveform. In Figure 2-15 the 
top trace shows the signal plus the burst noise, the bottom 
figure shows it after it has been maeenea at 100 hz. Fig- 
eres 2-16 and 2-1/7 show the burst noise alone in the top 
Berace, the signal plus noise Pacers: noise in the middle 
trace and the bottom trace is the middle trace after it has 
been filtered at 100 hz. In each figure the time scale is 


0.25 seconds per division. 
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As was previously doné the data was sampled at 256 sam- 

| ples per second and 2048 point transforms were then computed. 
Figures 2-185 through 2-21 show the results of the transforma- 
tions. Each trace contains the average of 40 transforms and 
as was done previously there is 80% redundancy between adja- 
cent traces. The S/N ratios on the figures refer to the S/N 
ratios of the signal to continuous noise and do not account 
for the burst noise. 

Figures 2-22 through 2-28 are photographs of the same 
material. The photographs are presented because the graph 
routine that drew Figures 2-18 through 2-21 would not func- 
tion if fewer than five traces were to be drawn. Figures 
2-22 through 2-24 present the transforms in the same manner 
as Figures 2-18 through 2-21. Each trace is the summation 
or average of 40 transforms. Figures 2-25 through 2-25 
present the results of averaging 80 transforms. Figure 2-25 
Shows the -24 db S/N transforms and Figure 2-26 is the same 
presentation only in this case the coefficients have been 
eaquared. It is felt that squaring a spectrum gives better 
insight into the Signal recognition problem. If the signal 
is readily discernable before the coefficients are Squared, 
then after squaring it will be even more discernable. If 
the signal becomes less discernable after squaring, then it 
is concluded that the S/N ratio is insufficient to allow any 
degree of Signal recognition. 

Figures 2-27 and 2-28 present the -30 db transforms 


exactly as described for the -24 db transforms. 
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Tables IV and V give the Signal to noise output ratios 
for the different S/N ratios. Table IV is the average S/N 
output ratio for the summation of 40 transforms, and Table V 
is for the summation of 80 transforms. In each of these 
cases the noise statistics were computed using 25 coeffi- 
cients on either side of the signal coefficient, not includ- 
ing the coefficients adjacent to the signal coefficient. 

The analysis of the meee Armat LORS presented in this 
section lead quite conclusively to the result that the modi- 
fied FFT is capable of providing better S/N output ratios 
than the FFT on this type of data. The modified FFT provides 
substantially improved S/N output ratios over the FFT in 
Figures 2-18, 2-19, and 2-20. In Figures 2-21 neither trans- 
formation provides sufficient signal recognition. 

Figures 2-22 through 2-24 are photographs of the spectra 
presented in Figures 2-18 through 2-20. The photographs 
emphasize the extra noise smoothing provided by the modified 
FFT. Figures 2-25 through 2-28 result from the average of 
80 transforms instead of 40 as was done in the previous fig- 
ures. In each of these figures the variance of the noise 
coefficients is considerably greater in the FFT spectrum 
than in the modified FFT spectrum. In each of the figures 
the signal is more distinct and more easily recognized in 
the modified FFT spectrum than in the FFT spectrum. 

Table IV shows the consistently better S/N output ratios 
of the modified FFT over the FFT. The -24 db input S/N 


ratio is the only exception. Another interesting observation 
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is the improvement of the S/N output ratios obtained by mul- 
eepiying thestr! and modified FFT spectra together. The 
product of the two spectra is consistently better than either 
the individual Spectra or the FFT spectrum squared. 

Table V further emphasizes the modified FFT's merits in 
transforming this data. By averaging 80 transforms the mod- 
ified FFT was able to produce better S/N output ratios than 


the FFT in every case tested. 
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40 transforms averaged per trace, transform length 208 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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hO transforms averaged per trace, transform length 208 
Displayed bandwidth 18.75 hz, 0.725 hz resolution 
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4O transforms averaged per trace, transform length 20)8 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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hO transforms averaged per trace, transform length 208 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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Figure 2-19 (Mod-FFT) 
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hO transforms averaged per trace, transform length 20,8 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


ceweome = O06 oi WNT odie, 
feleeee = J.99—8-41 UNITO7 INCH. 


S4N45URS7T FRE 
BY N=-04 DR Figure 2-20 (FFT) 


83 











hO transforms averaged per trace, transform length 20)8 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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Figure 2-20 (Mod-FF™) 
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hO transforms averaged per trace, transform length 20)8 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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4O transforms averaged per trace, transform length 20)8 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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Figure 2-21 (Mod-FFT) 
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Modified FFT 
10 hz signal, S/N = -12 db + burst noise 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-22 
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Modified FFT 


10 hz signal, S/N = -18 db + burst noise 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 
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Modified FFT 


10 hz signal, S/N = -24 db + burst noise 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-24 
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Modified FFT 


10 hz signal, S/N = -24 db + burst noise 
80 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


Figure 2-25 
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Modified FFT 
10 hz signal, S/N = -24 db + burst noise 
80 transforms averaged per trace, transform length 2048 


Displayed bandwidth 18.75 hz, 0.125 hz resolution 
Spectrum squared 


Figure 2-26 
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Modified FFT 
10 he signal, S/N = -30 db + burst noise 
80 transforms averaged per trace, transform length 2048 
Displayed bandwidth 18.75 hz, 0.125 hz resolution 


| Figure 2-27 
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Modified FFT 


10 hz signal, S/N = -30 db + burst noise 
80 transforms averaged per trace, transform length 2048 


Displayed bandwidth We. 75° hz, 0. L2oe nz resolution 
Spectrum squared 


Figure 2-28 
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S/N out (db) 


a. FFT Mod-FFT 9 FPTXFFT 4 ipa 
-12 13.770 22.354 16.342 23.711 
-18 7.729 14.389 8.671 15.133 
-24 8.747 8.360 6.450 8.863 
-30 5.109 5.515 3.551 5.910 

Table IV 


Average S/N out (db) vs. S/N in (db) for the Average of 40 
Transforms for the case of the 10 hz signal plus Gaussian 
noise plus burst noise. 
_ 2 
(Fo-u,) 


PX oi 
n 


S/N out computed as: 


where Fe is the magnitude of the signal coefficient and HL 
and Oo, are computed from 25 noise coefficients on either 


Side of F._. 
S 
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oy Neen 
(db) 


~12 
-18 
-24 


-30 


Average S/N out (db) vs. 


(ak 


16.154 


11.433 


9.332 


6.875 


S/N out (db) 


Mod-FFT 


26.606 


16.341 


m0 2 19 


8.094 


Table V 


dod hese daly 


18.731 


12.855 


6.713 


9.130 


FFTx 
Mod-FFT 


26.945 


lg salpa 


10.061 


(eccloab 


S/N in (db) for the Average of 80 


Transforms for the case of the 10 hz signal plus Gaussian 


noise plus burst noise. 


Cte 


D6" 
n 


S/N out computed as: 


where F. is the magnitude of the signal coefficient and WD 


and Oo, are computed from 25 noise coefficients on either 


side of F.. 
S 
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E. UNDERWATER ACOUSTIC DATA 

The final data to be analysed was underwater acoustic 
data. This data provides a great variety of signals from 
different sources with different strengths all transmitted 
Simultaneously. It is a challenging test of the signal proc- 
essor's ability to resolve the many different signals and to 
perceive weak signals in the presence of strong ones. Only 
selected portions of the data have been presented in this 
section. 

The data was filtered at 200 hz and sampled at 512 sam- 
ples per second. Again, 2048 point transforms were employed 
resulting in a frequency resolution of $ hz. 

Each trace in the pictures is the result of averaging 
40 transforms. The FFT spectra are in the top photographs 
in each figure and the modified FFT spectra are in the lower 
photographs. The time frame and frequency bandwidth are the 
same for each photograph in a given figure. The total dis- 
played frequency bandwidth in each trace is 37.5 hz. 

Figures 2-29 and 2-30 each depict the same frequency 
bandwidth but at different times. The main distinction is 
the modified FFT's ability to detect the highlighted signal 
in Figure 2-29 before the FFT and to make it clearly recog- 
nizable. The large signals to the left and right of the 
highlighted signal are enhanced better by the FFT. There 
are also significant differences in the general appearance 
of the two spectra. Figure 2-30 is presented to show that 


at a later time the FFT presented the same signal as well as 
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or better than the modified FFT. The signal to the left of 
the highlighted signal seems to be presented equally well by 
both spectra. 

Figures 2-31 and 2-32 show a relatively weak signal. 
Close comparison of the traces in the spectra indicate that 
in some cases the FFT outperforms the modified FFT and in 
others the reverse is true. The bottom trace of Figure 2-32 
provides a dramatic case in which the modified FFT produces 
the signal better than the FFT. 

Figures 2-33 through 2-35 provide an example in which 
the two spectra behave almost identically except for the 
last trace in Figure 2-35 in which an effect like the one 
described for Figure 2-32 above occurs. 

Figure 2-36 again depicts a Signal that has been en- 
hanced better by the modified FFT than by the FFT. Also two 
distinct signals are present to the right of the highlighted 
Signal in the modified FFT whereas only one is present in 
the FFT. 

Figures 2-37 through 2-39 all show the same signal at 
different times. In this case the two spectra are again 
quite Similar. The main reason for presenting these spectra 
is the modified FFT's ability to retain the signal in Figure 


2-39 after it has been lost by the FFT. 
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Modified FFT 
Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-29 








Modified FFT 


Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-30 
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Modified FFT 


Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-31 
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Modified FFT 


Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-32 
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Modified FFT 
Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-33 
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Modified FFT 
Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-34 
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Modified FFT 
Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-35 
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Modified FFT 


Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-36 
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Modified FFT 


Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


FIrgquce 2-37 
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Modified FFT 


Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-38 
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Modified FFT 
Underwater Accoustic Data 
40 transforms averaged per trace, transform length 2048 
Displayed bandwidth 37.5 hz, 0.25 hz resolution 


Figure 2-39 
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IV. PART 3 SUMMARY OF EXPERIMENTAL RESULTS 


The data analysed in this paper were basically of three 
types. A signal or signals with additive Gaussian noise, a 
Signal with additive Gaussian noise plus bursts of Gaussian 
noise spaced at random intervals with random duration, and 
finally some underwater acoustic data. 

In Section III-B of this investigation, a single sinu- 
soid with additive Gaussian noise was processed with the FFT 
and the modified FFT. Figures 2-1 through 2-9 indicated 
that the resulting spectrums from the two transformations 
are quite Similar. The FFT seems to have smoothed the noise 
slightly better than the modified FFT. In general the FFT 
seems to provide the better processing gain, but upon close 
scrutiny of the spectra it iS possible to find some traces 
in which the modified FFT seems to produce a greater signal 
to noise ratio. The computed S/N output ratios in Table I 
seem to agree quite well with what is observed in the photo- 
graphs. The modified FFT seems to outperform the FFT at the 
large S/N ratios and the FFT dominates the modified FFT at 
the small signal to noise ratios. 

Section III-C of this paper dealt with multiple signal 
data with additive Gaussian noise. Figures 2-10 through 
2-14 depict the output spectra and again it is noted that 
the spectra are quite Similar in overall appearance, with 
neither transformation appearing to have a distinct advan- 


tage over the other. The most interesting results are in 
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Figures 2-13, where in the FFT spectra the -16 db signal is 
consistently smaller than the -14 db signal. In the modi- 
fied FFT spectra the -16 db signal is nearly as large as the 
-14 db signal. From the -12 db signal on down it is doubt- 
ful that one could recognize the signal presented in the 
figures. Longer transforms or the averaging of more trans- 
forms would be required to enhance the signal to noise ra- 
tios. Averaging twice as many transforms theoretically 
increases the S/N output by 3 db and this is verified in 
Tables II and III. 

These tables tend to indicate that the FFT is the better 
transformation at high S/N ratios and the modified FFT is 
better at low S/N ratios. For reasons pointed out in Section 
IITI-C the -18 db Signal should not be used in the analysis. 

Some of the apparent inconsistencies in thesSe tables are 
explained by the fact that different noise coefficients were 
used in computing the noise statistics for each Signal as 
mentioned in Section III-C. There is no known explanation 
for the inconsistency of the process gain for the -6 db Sig- 
nal. A search through the adjacent coefficients showed that 
this coefficient provided the best signal to noise output in 
the vicinity of the signal for the FFT. 

Another interesting aspect of Tables II and III is that 
in going from a 40 transform average to an 80 transform aver- 
age the modified FFT S/N ratios are improved more than the 


FFT S/N ratios. 
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section III-D presents data in which bursts of noise 
were added to the signal in addition to continuous Gaussian 
noise. In this case the modified FFT consistently outper- 
formed the FFT. The one exception was in the case of the 
-24 db signal where 40 transforms were averaged. When 80 
transforms were averaged the modified FFT outperformed the 
FFT in every case. Figures 2-25 through 2-28 show the re- 
sults of the 80 transform average in which it is much easier 
to believe in the existence of a signal in the case of the 
modified FFT than in the FFT. 

Tables IV and V present the S/N ratios for the spectra. 
Little needs to be said in this case since the results seem 
to be consistent with the observations made in Figures 2-18 
through 2-28. 

section III-E of this paper presents Several Signals en- 
countered in the transformation of some actual underwater 
acoustic data, an area in which signal processing finds a 
very practical use. The manner of interpretation of the data 
is very subjective, as is the manner in which it is normally 
analyzed. In the signals that were presented the modified 
FFT usually did as well or better than the FFT in making sig- 
nal presence easily detectable. During all of the research 
that was done it was found that both transformations gener- 
ally provided similar results for this type of data. Ex- 
tremely strong signals were generally enhanced more by the 
FFT. Medium to weak signals appeared to be about the same 


in overall detectability for both transformations. Looking 
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specifically from trace to trace differences were definitely 
noted and no one particular transformation proved to be con- 
sistently better. On some occasions the modified FFT pro- 
duced signals that the FFT could not, even eee later on 
the signal was produced by the FFT. In the case of Figure 
2-29 an analyst confirmed the high probability of a signal 
being present at the spotlighted frequency. In Figure 2-39 
the ability of the modified FFT to retain the signal longer 
than the FFT is another strong point in favor of the modi- 
fied FFT but an unusual occurrence. The nonlinearity of the 
modified FFT with the transformation being a function of the 
input data itself has made it difficult to draw generalities 
about its performance except possibly in the case of data 


Similar to that presented in Section III-D of this paper. 
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V. PART 4 CONCLUSIONS 


The data analyzed in this paper is far from conclusive 
or absolute in terms of allowing any hard and fast conclu- 
Sions to be drawn. The process does seem to have some mer- 
its, especially in the case of the signal plus noise plus 
burst noise. Whereas this may not appear to be a very prac- 
tical type of data, it is felt that somewhere between this 
case and the case of stationary Gaussian noise a practical 
Situation exists in which this transformation may be very 
useful. Somewhere where strong transients exist, whether 
they are switching noises, or the explosion of underwater 
sonic sources, this transformation may be useful if an under- 
lying continuous signal is of interest. 

Unexplained are the instances in Section III-E of this 
paper where the modified FFT was able to bring out Signals 
Significantly better than the FFT, or in one case to retain 
the signal for a longer period of time than the FFT. 

On the other hand, the FFT seems better able to distin- 
guish between signal and noise when the noise is stationary 
Gaussian. 

The good characteristics of the transformation do not 
seem to be dependent on the signal frequency. In Section 
III-E the signals that were much more evident in the modi- 
fied FFT than in the FFT, appeared at both ends of the 
spectrum as well as in the middle of the spectrum. This ef- 
fect was not limited to one data tape but appeared on two 


different tapes from different sources. 
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The main drawback to implementing the modified FFT as it 
is presently being utilized is the length of time required 
to perform the transformations, approximately three to four 
times as long as a similar FFT transformation for 1024 points. 
However if parallel processing were available the modified 
FFT could be performed twice as fast as the FFT since each 
pass only operates on half of the coefficients. 

In general there does not — to be enough difference 
between the FFT and modified FFT spectra to warrant its use 
in preference to the FFT. If it is to be used it would 
seem appropriate to use it when a signal is suspected and 
the FFT does not reveal it. Also, if the data were known to 
have characteristics similar to that of the signal plus noise 
plus burst noise it would warrant strong consideration. 

One other asset of the transformation that may prove impor- 
tant is the inability of any given output coefficient to 
exceed some maximum value, determined by the transform length. 

With so many variables, such as transform length, fre- 
quency resolution, windowing, Signal stability, and differ- 
ing types of background noise and signals, much more 
investigation will be required to determine the optimum 
utilization of this transformation. A mathematical and 
Statistical analysis of the modified FFI may give further 


imc tS best applicavion. 
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APPENDIX A 


FAST FOURIER TRANSFORM - 


The Fast Fourier Transform is an algorithm for computing 
the "Discrete Fourier Transform" which evaluates the Fourier 
Coefficients of the Fourier Series for a discretely sampled 
waveform. 

The term Fast Fourier Transform is somewhat misleading 
because the FFT is applied to the transformation of periodic 


functions which do not necessarily satisfy the condition: 


co 


pret) | enero 


which is required for f(t) (in general) if it is to have a 
Fourier Transform. 

Any discrete analysis of a function f(t) must of neces- 
sity be finite in length, say from t=0 to t=T. What happens 
before t=O and after t=T is in general not known, or at 
least the time record in question gives no indication of 
what was occurring prior to t=0 or what is going to occur 
after t=T, unless the function iS periodic and the time re- 
cord contains at least one whole period of the waveform. 
When dealing With the FFT the implied assumption is that the 
record to be analyzed is one period of a periodic function 
and that the information contained in this period completely 
describes the function from t=-© to t=to, 

Periodic functions are more appropriately described by 


PiemwOourier oeries expansion, rather than by the Fourier 
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Transform. The term Fast Fourier Transform may be more de- 
scriptive if it were renamed the Fast Fourier Series. 

A waveform f(t), periodic in t with period T satisfying 
£(t)=f£(t+T) for all t with a finite number of discontinuities 
on the interval (0<t<T) can be expanded in a Fourier Series 
as: 

f(t) = = F(n)exp[((j2mnt)/T} 

n=— 
where F(n), the Fourier Coefficients are defined as: 
i | 

F(n) = 7 : f(t)exp[(-j2nnt)/T]dt 

The above representation of the Fourier Coefficients 1s 
for the continuous case. In order to determine the coeffi- 
cients for the discrete case the integral must be approxi- 
mated by a summation. Assuming the waveform has been sampled 
N times in T seconds and that the samples are spaced AT sec- 
onds apart, then, 

T=NAT 
and 

t=kAT (k=0,1,2,...,N-1) 


then the sampled version of the Fourier Coefficients can be 


written: 
1 N-1 
F(n) = yam © f(kAT)exp((-j2akAT)/T]AT 
k=0 
1 N-1 
=— 2 f£¢kAT)exp[(-Jj2mnkT)/ (NT) ] 
k=0 
1 N-1 
=5 2 £(kAT)exp[(-j2mnk)/N] (n=0,1,2,...,N-1) 
k=O 
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The factor 1/N is a scaling factor, useful for obtaining 
the original amplitude of the Apa waveform when f(t) is 
expressed in its Fourier Series. It has no effect on the 
relative amplitudes of the Fourier Coefficients. 

F(n) is the amount of frequency n/T contained in f(t). 
Abbreviating f(kAT) as f(k), and omitting the scaling factor 


1/N for convenience, then, 


N-1 
F(n) = ££ f£f(k)exp[(-j2mnk)/N] 

k=0 
N-1 

= Y f(k){cos[{(2tnk)/N]-jsin[{(2mk)/N]} 
k=0 
N-1 N-1 

= 2 f(k)cos[((2mnk)/NjJ-j «© f(k)sin((2mnk)/N}] 
k=0 k=0 


which is simply a discrete correlation of f(t) with N sinu- 


soids at the frequencies n/T Since, 


N NAT ~ Tar ~ ae 
F(O) is then the amount of D. C. component contained in f(t) 
F(1), the amount of frequency 1/T = Af contained in f(t), 
F(2), the amount of frequency 2/T = 2Af contained in f(t), 
etc., up to n=N-1 for a complex waveform f(t). 

If f(t) is real valued the Fourier Coefficient F(N/2-m) 
= F’(N/2+m) (where * denotes the complex conjugate). 


Proof: 
N-1 

y f(k){cos[(2mnk)/N)-jsin[ (21nk)/N] } 
k=0 


It 


Z| 


F(n) 


jet Fin) = R(n)-jit€n), then, 
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N-1 


R(n) = Re[F(n)] = = & £(k)cos[(2mnk)/N] 
k=0 
| Ned ie 
for n=N/2, R(N/2)= N x f(k)cos(mk)= F y (-1) £(k) 
=0 k=0 


now: cos[27Tk(N/2-m)/N] = cos[{(Nrk-27Tkm)/N] 


fl 


cos[(Ntk)/N]cos[(21Tkm)/N] + sin[(Nrk)/N]sin[(27km)/N] 
where the second term is zero for integer k. 

= cos[(Nkt + 21km)/N] = cos[27k(N/2 + m)/N] 
therefore, 

R(N/2-m) = R(n/2tm) (m=0,1,2,...,N/2-1) 
or the real part of F(n) is symmetric about n=N/2. 

A similar proof shows that the imaginary part of F(n) is 
antisymmetric about n=N/2, i.e., 


I(N/2-m) = -I(N/2tm) (m=0,1,2,...,N/2-1) 





Figure A1l-1 


That is F(N/2-m) = F'(N/2+m), mpivane chat all of the 
spectral information is contained in the first N/2 coeffi- 
cients, the second half of the spectrum containing redundant 
information. This is intuitively obvious if one considers 
that the cosine is an even function and the Sine is an odd 


huner lon. 


118 





To more readily visualize the correspondence between the 
Fourier Transform of a continuous and a sampled waveform, 
consider a continuous waveform f(t), and its samples counter- 
part T(t) (where the k denotes the sampled version of f(t)). 


Then f(t) can be represented as: 


f(t) = f¢(k) = %& f(t)dé6(t-kAT) 
K=-0 
let d(t) = %& 6(t-KAT) 
k=-—c 
Diw) = d(t) 
F(w) = f(t) 


FCW) <> f(t) 


where implies Fourier Transform pairs. 


Then FW) = D(w)*F(w) 


Since {6(t-kAT)} is a train of impulse functions periodic 
in t with period AT, the Fourier Series representation is 


used to describe this function, that is: 


co 


d(t) = %Y D(n)exp(jwnt) where w=2nf = 21n/AT 
n=-© 
and 
1 AT 
D(n) = AT f d(t)exp(-jwnt)dt 
0 
then, 
1 AT 
Din) =— f «: 6(t-KAT)exp(-jwnt )dt 
i 
a °° 
=e Oy f 6(t-kAT)exp(-jwnt)dt 
Lt ao! 





oe Le. ¢) 


1 : 
D(n) = AT © exp(-jwntkAT) = mei ee xD Jeera | 
k=— 0 k=-0 


jas 4 
aT,” exp(-j2mnk) = AT for all integer n, k 


— © 


Now, 


6(t-kAT) = —= © exp(jwnt) 


n=-© 


d(t) 


then: 


© 


s{f(t oe » exp(jwnt ) jexp(-jwt)dt 


n=—0o 


FL Cw) 


=s 5 Sf{f(t)exp (jwnt) Jexp(-jwt)at 


N=-©0 -—© 


© foe] 


a Sf £(t)exp[jCwn-w)t jdt 


n=--o —-© 


i 
AT 


© © 


» Sf £(t)exp[ 3 (Cwn-w)t jdt 


n=-© -—-o 


_i 
AT 


= — ¥ F(Cw-wn) 


The discrete Fourier Transform: 
F(n) = = 5 £(kAT)exp(-j2mnk/N) (n=0,1,2,...,N-1) 


can be visualized in various ways. AS previously mentioned, 
it can be considered a discrete correlation of the sampled 
waveform f(kAT) with a Series of sSinusoids at the frequencies 
Mieco Oem ere N-1) where T=NAT is the length of the input 
waveform to be transformed, N is the total number of samples, 


and AT is the time between samples. 


1.20 








The time T then determines the frequency resolution of 
the transformation, that is, 1/T=Af or the difference in 
frequency from one correlating Sinusoid to the next. N and 
T determine the highest frequency of the correlating sinu- 
soids. The ratio of N/T is the sampling rate, i.e., N sam- 
ples in T seconds. For real valued input functions, the 
highest correlating frequency is then (N/2-1)/T or half of 
the sampling rate. Figure Al-2 shows that the waveform must 
be sampled at a rate at least twice as high as the highest 
frequency present in f(t) or the spectra of F(w) will over- 
lap. This is called aliasing. 

A second visualization is simply that of a matrix opera- 
tion. Expressing the discrete Fourier Transform in matrix 
form, omitting the scaling term 1/N, and defining W as exp 


(-j27/N) then: 


[F(n)} = (WwW) £Ck)] 
An expansion of these matrices for eight points is given in 
Figure Al-3a. 
Figure A1l-3b shows the matrix rewritten with the follow- 
ing substitutions: 


nk nkModN 


W = W Since W is periodic with period N. 
Also, W° = exp(0) = 1 = -W"* 
W? = exp(-j7r/4) = -W° 
W2 = exp(-j1/2) = -W® 
W? = exp(-j37/4) = -W’ 


The matrix in Figure A1-3b can be expanded into three 


matrices as shown in Figure Ai-4. 
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Resulting spectrum if f(t) has been sampled 
at a rate at least twice as fast as the high- 


est frequency present in f(t). 
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Resulting spectrum if f(t) has been sampled 


at a rate slower than twice the highest 


frequency present in f(t) 


Figure Ai-2 
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Discrete Fourier Transform in matrix form 


for 8 points 
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Matrix in (2) rewritten usin the substitution 


yok = 


nkModN 


Figure At-3 
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It is this matrix expansion, along with the periodic 
properties of yak (the correlating sinusoids), that gives 
rise to the Fast Fourier Transform. In general, if a sam- 
pled input waveform has N samples where n=2™ and M is a pos- 
itive integer, then the NXN matrix of operators (the matrix 
of W's) can be expanded into a product of M matrices with 
the stipulation that the input vector is bit reversed (note 
the reordering of the input vector). This stipulation is 
not a requirement but it is the method that is used through- 
out this study. The process of bit reversal allows the ma- 
trix multiplications to be performed one after the other 
without a requirement to save any intermediate results. The 
matrix expansion can be performed in other ways such that 
the input vector is left in its natural order and the output 
vector is in bit reversed order, or possibly the odd harmon- 
ics of Af are grouped together and the even harmonics are 
grouped together. It is also possible to perform the opera- 
tions without bit reversing the input or the output vectors 
but then a requirement exists for the retention of inter- 
mediate results. 

Bit reversal is simply the reversing of the order of the 
digits that make up a number in binary arithmetic. In the 
above situation N=8, and M=3. For example, 110 (6) becomes 
011 (3) when its bits are reversed. The effect of bit re- 
versal is to place each element in the input vector next to 
an element that is N/2 elements away from it, i.e., N/2=4 in 
this case so element 1 is placed next to element 5, 2 next 


tomo, etc. 
ZS 





To go one step further, the matrices of Figure Al-4 can 
be further reduced by using a tree graph or what is commonly 
called an "FFT butterfly". The expanded matrices have many 
zero elements, in fact each matrix contains only 2N elements 
that are non-zero and a large portion of these are 1's or 
-1l's. The original Discrete Fourier Transform required N? 
complex operations in its execution; the expanded matrices 
of the Fast Fourier Transform only have Nlog.N operations 
and thus the term ''Fast". 

The matrices have been rewritten in butterfly format in 
Figure Al-5. The column labeled "0" is the input vector. 
The column labeled "1" is the bit reversed form of the input 
vector. The column labeled "2" is the results of the first 
matrix multiplication, or the results after the first pass 
through the algorithm. The number 2 is used to describe 
this column since two of the original elements make up each 
of the elements at this point in the algorithm. The column 
labeled "4'' is the result of the second matrix multiplica- 
tion, or the result after the second pass through the algo- 
rithm. As above the number 4 is used to describe this 
column since four of the original elements make up each of 
the elements at this point in the algorithm. The nomencla- 
ture is similar for each of the columns in the butterfly. 

When two arrows meet at a circle with a + or - in it, 
then addition or subtraction is to take place. For example, 
in the column labeled 2 element f) is made up of elements 


£5 + f£, from column 1; element f; in column 2 is made up of 
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FFT butterfly or flow diagram for an 8 point 


transformation 


Figure A1-5 
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elements fo - f£,; from column 1. The circles that contain an 
argument like 1, 2, or 3 imply a multiplication of the pre- 
ceding element by W’, W’, or fEmbocoroeen eddition or Sub— 
traction is performed. For example, in column 4 element 

fo is made up of elements fo + f2 from column 2, element 

f; is made up of elements f; + f3W* from column 2, element 
fe is made up of elements fo - f2 from column 2, element 

f3; 1s made up of elements f; - f3W* from column 2. 

The butterfly is identical to the matrix approach in 
that each pass through the butterfly is identical to perform- 
ing one of the matrix multiplications in the matrix repre- 
sentation of the FFT. The butterfly condenses the matrix 
approach by elimination of the zeros in each matrix leaving 


only the necessary operations. 
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